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ABSTRACT 

We present a new analysis of stellar mass functions in the COSMOS field to fainter limits than has been previ- 
ously probed at z < L The increase in dynamic range reveals features in the shape of the stellar mass function 
that deviate from a single Schechter function. Neither the total nor the red (passive) or blue (star-forming) 
galaxy stellar mass functions can be well fit with a single Schechter function once the mass completeness limit 
of the sample probes below ~ 3 x lO^M©. We observe a dip or plateau at masses lO'^M©, just below the 
traditional M* , and an upturn towards a steep faint-end slope of a ~ -1.7 at lower mass at all redshifts < L 
This bimodal nature of the mass function is not solely a result of the blue/red dichotomy. Indeed, the blue 
mass function is by itself bimodal at z L This suggests a new dichotomy in galaxy formation that predates 
the appearance of the red sequence. We propose two interpretations for this bimodal distribution. If the gas 
fraction increases towards lower mass, galaxies with Mbaryon ~ 1O'"M0 would shift to lower stellar masses, 
creating the observed dip. This would indicate a change in star formation efficiency, perhaps linked to super- 
novae feedback becoming much more efficient below ^lO'^M©. Therefore, we investigate whether the dip is 
present in the baryonic (starsH-gas) mass function. Alternatively, the dip could be created by an enhancement 
of the galaxy assembly rate at ~1O"M0, a phenomenon that naturally arises if the baryon fraction peaks at 
Mhaio IO'^Mq. In this scenario, galaxies occupying the bump around would be identified with central 
galaxies and the second fainter component of the mass function having a steep faint-end slope with satellite 
galaxies. The low-mass end of the blue and total mass functions exhibit a steeper slope than has been detected 
in previous work that may increasingly approach the halo mass function value of -2. While the dip feature is 
apparent in the total mass function at all redshifts, it appears to shift from the blue to red population, likely as a 
result of transforming high-mass blue galaxies into red ones. At the same time, we detect a drastic upturn in the 
number of low-mass red galaxies. Their increase with time seems to reflect a decrease in the number of blue 
systems and so we tentatively associate them with satellite dwarf (spheroidal) galaxies that have undergone 
quenching due to environmental processes. 

Subject headings: surveys — cosmology: observations — galaxies: mass function — galaxies: evolution 



1. INTRODUCTION 
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Galaxy formation and evolution is believed to be driven pri- 
marily by two processes: firstly, the successive merging of 
their parent dark matter halos causing accretion of material 
and ultimately mergers between galaxies; and secondly, the 
feedback-regulated conversion of gas into stars within galac- 
tic disks with subsequent potential rearrangement of the disk 
material by dynamical processes (secular evolution). Both 
processes contribute to the growth in stellar mass of galaxies 
with time. The stellar mass function of galaxies and its evolu- 
tion with time is therefore fundamental to the understanding 
of galaxy formation. 

The ability to estimate galaxy stellar masses has ad- 
vanced in recent years in large part because of increas- 
ing access to near-IR photometry. Estimates are typically 
made by fitting multi-band photometry with stellar pop- 
ulatio n synthesis librarie s (see , e.g., Brinchmann & ElliJ 
1200^, 'Bruzual & Chariot' "2003^, "Drorv et al. 2003, 2004 
Maraston et al, 2006; Marchesini et al. 2008; Conrov et al] 
2009), fitt ing specific spectral fea tures when spectroscopy is 
available jKauffmann et al.|[2003h . or the full spectrum when 
high-quahty spect ra are observed (iReichardt et al.l 120011: 
iPanter et"aLl 120041) . So far, only the photometric fitting 
technique has been a viable option for high-redshift sur- 
veys. These measurements provide masses with accuracies 
of ~ 0.1-0.3 dex, and systematic uncertainties of up to 
a factor of two depending on the selection and number of 
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photometric bands included and the assumptions made on, 
among others, the shape of the IMF, the allowed star for- 
mation histories, the dust extinction model, or the underly- 
ing st e llar population synthesis method (see, e.g.,[P rorv et al. 
2004; Kannappan & Gawisei' '2007, 'M archesini et al/BOOSl: 
Conroy et al. 2009 for systematic studies on this matter). 

Utilizing such techniques, the build-up of the stellar mass 
density from redshift z ^ 6 to the present epoch has been the 
subject of several studies in the past decade, often relying on 
deep multi-band imaging surveys in the UV to near-infrared 
wavelength range (Brinchmann & Ellis 2000; Drorv et al. 
2001; Cohen 2002; Dickinson et al. 2003; Fontanaetal. 
2003; Rudnicketal. 2003; Glazebrook et al. 2004; 
Cpnselice et al. 2005; Chapman et al. 2005; Dro rv et alJ 
2005(; Rudnick et al. ^006t lEvles et al.i i2007c iGrazian et alj 
200l [stark et al. 200l. 



The stellar mass function in the local universe has been 
measured from large imaging and spectroscopic survey s such 
as 2dF, SPSS, and 2M ASS (Cole e t ani200l1; Bell et al .ll2003t 
Perez-Gonzalez et al.i 12003; .PanterTt all 12004; .Baldrv et alJ 
20081) . At distances up to z ~ 1.5, a number of groups 



have established a picture of the evolution of the mass 
function with some detail (Drorv etal. 2003; Fontana et aL 
■2004; Bundvetal. 2005; Borchetal. 2006; Bund v et al 
2006; Arnouts et al.. .2007; Pozzetti et al.. .2007; Ilbert et al 



2009bl) . with generally good agreement between different 



datasets. To some lesser detail and accuracy, deep surveys 
have provided data spanning < z < 5 (Drorv et al. 2005; 



Cpnselice et al." '2005; Fontana et al. 2006; Yan et al. 2006; 
Grazian et al. 2007; Eisner et al. 2008; Perez-Gonzalez et al, 
2008r Marchesini et al. 2008), and even some estimates at 



z ~ 7 (iBouwens et al..i200 6). So far, these high-redshift stud- 
ies of the stellar mass function emphasized the evolution of 
galaxies of mass logM/M© > 10. Speaking very broadly, the 
stellar mass density decreases by a factor of two to z ^ 1, with 
the most massive galaxies already being in place at earlier 
epochs. The evolution appears to accelerate quickly beyond 
z- 1.5. 

In this paper, we concentrate on the low-mass galaxies that 
have typically been below the completeness limits of pre- 
vious work at z > 0.1. Generally iSchechterl (11976) fits to 
the galaxy stellar mass function with faint-end slope -1.1 
to ~ -1.3 have been found adequate to describe the galaxy 
population (even separated morphologically, by color, or star 
formation activity; Pannella et al. 2006; B orch et al.l 120061; 
lArnouts et al.1l2007t IPannella et al.ii2009i: illbert et al.ll2009bh . 



Recently, though, a steepening of the slope of the lumi- 
nosity function below M, ^ -17 in the local universe has 
been convi ncingly detec t ed in clusters ([Driver etaTI 119941; 
Trentham & Tullvl 120021: iHilker et al.1 12003 1 ; iPopesso etakl 
2005, 2006), groups (Trentham & Tully"2002';'Trentham et al.) 
2005; Gonzalez et al. 2006), and in the field (Blanton et al.' 
20051) . For example, iTrentham & TuUyl (120021) find that the 



luminosity function in the Virgo cluster, in the NGC 1407 
group, and in the Coma 1 group is steep between Mr of - 1 8 
and -15 (and flatten s again only at Mr > -15). Moreover, 
iBaldrv et al.l (120081) find that the local galaxy stellar mass 
function steepens as well below logM/M© ~ 9.5 (but see also 
ID & White 2009). 

The steepening of the mass function can also be interpreted 
as a bimodality: the mass function consists of a sum of (at 
least) two components. This bimo dal behavior has no w also 
been detected at redshifts z > 0.1. iPozzetti et aTl ('2009) find 
bimodal mass functions to z ~ 0.5 from the zCOSMOS spec- 



troscopic survey. They interpret the mass function as being 
composed of early-type galaxies dominating the massive part 
and late-type galaxies dominating the less massive part and 
contributing the steep faint-end slope. Each of these com- 
ponents is well fit by a Schechter function. Bolzonell a et aP 
(12009) use the same sample to investigate the bimodality as a 
function of environment. They find that at z < 0.5, the shape 
of the galaxy stellar mass function in high- and low-density 
environments become markedly different, with high density 
regions showing a stronger bimodality. 

We extend the study of the shape of the galaxy stellar mass 
function, particularly at low masses, to z 1, with stellar mass 
limits ^1.5 dex lower than can be achieved with spectro- 
scopic studies. We investigate actively star-forming galaxies 
and passively evolving galaxies separately; we study how the 
change in slope may be caused by the presence of multiple 
galaxy populations, that taken together, lead to a mass func- 
tion shape that is more complex than a single power law with 
an exponential cutoff or even a simple combination of early- 
and late-type components. We show that the blue mass func- 
tion itself is bimodal and that passive galaxies exhibit a faint- 
end upturn, likely caused by dwarf spheroidal galaxies linked 
to the faint end of the blue galaxy population. 

This paper is organized as follows: in §|2]we introduce the 
galaxy sample that we use in this work. In §[3] we discuss the 
stellar population models used to derive stellar masses and 
the resulting mass completeness limits. In §|4]we present the 
stellar mass function of active and passive galaxies and we 
discuss our results in § |5] Finally, we summarize this work in 

m 

Throughout this work we assume J7m = 0.3, J7a = 0.7, //q = 
70 /lygkm s"' Mpc"'. Magnitudes are in the AB system. We 
will denote galaxy stellar masses by the symbol M - or Mg 
where an explicit distinction form halo masses, denoted by 
Mh, is necessary. The symbol M* is reserved for the charac- 
teristic mass parameter of the Schechter function. 

2. SAMPLE SELECTION 

The primary focus of this paper is determining the abun- 
dance of low-mass galaxies and characterizing the shape of 
the mass function, especially at the low-mass end. These 
galaxies are faint by definition, with / > 23 and K > 23, 
sources that are typically beyond the range of magnitude- 
limited spe ctroscop i c surv e ys like DEEP2 , WD S , or zCOS- 
MOS (Dav is etalJ 120031; iLe Fevre et all l200l iLiUv et al.1 
|2007|) . As such, we must rely on photometric redshifts. 

We use the COSMOS catalog with photometric redshifts 
derived from 30 b road and medium bands described in 
Illbert et all (l2009al) and ICapak et al.1 ( 120071) . version 1.5, 
dated April 05, 2008. We restrict ourselves to objects with 
i^g < 25. At this limit, the rms photo-z accuracy at z < 1 is 
~ 0.03 in Az/(1 +z) (Fig. 9 in Ilbert et al. 2009a). At z > 1 
the quality of the photometric redshifts quickly deteriorates. 
The d etection completeness at /^^ = 25 is > 90% ( C apak et al] 
I2007h . Additionally, we require a detection in the Ks band 
to ensure that the stellar mass estimates at the faint end are 
still reliable, thus we add the constraint Ks < 24. The surface 
brightness sensitivity limit for the /"^ data is 28.4 mag arcsec"' 
at 5(7. At this limit, an object of 25th magnitude would have 
to be 5.4 arcsec across (with no bulge) to have a chance of 
being missed. At 23rd magnitude that number increases to 7 
arcsec. It is hence very unlikely that surface brightness is the 
dominant selection effect at low luminosities where galaxies 
are very compact (in contrast to the situation in local-universe 




Fig. 1 . — A comparison of the selection techniques of passive galaxies in the NUV-r vs. r-J color-color plane. Galaxies identified as passive by SED fitting 
are marked as red points. The lines show the color cuts defined by WilUams et al. ( 2009) for identifying passive systems. In addition, galaxies with masses 
between the mass limit and 3 times the mass limit in each redshift bin are marked in blue, to look for low-mass galaxies with unusual colors which would indicate 
a possible problem with the photometric redshifts. 



surveys; see simulations in lCapak et al.ll20()7h . 

To ensure high-quality photometric redshifts we will limit 
the redshift range to z < 1 and to sample enough volume to 
z > 0.2. We have made sure that our results are not sensi- 
tive to this selection through repetition of our analysis with 
relaxed magnitude limits. Our mass functions do not change 
shape, nor change their faint-end slope if we drop the A'i-band 
selection and push the t selection to 26th magnitude. Drop- 
ping the Ks-hand selection leads to a gain in about one bin 
in depth (0.25 dex in mass), albeit with increased stellar mass 
errors by about 20% on individual sources below the Ks-hand 
limit, but no change in the faint-end slope (the z and /-band 
still provide enough information on the NIR light). Nor do 
our results change if we impose a brighter selection cut of 
i'^g < 24, apart from a restriction in redshift range. 

We use the photometric data in the u* (CFHT), Bj, Vj, g^, 
r+, /+, z+ (Subaru), J (UKIRT), and Ks (CFHT) bands for our 
analysis, totaling 9 bands. We explicitly exclude the IRACl- 
channel data (3.6 /im) due to their confusion-limited nature 
at niAB = 24 which would otherwise restrict the analysis to 
brighter magnitudes than is possible with the jT^ band data. 
This is i n contrast to the mass function analysis by lllbert et alJ 
(I2009hh . who use an IRACl -selected catalog. However, we 
have verified that adding the IRACl data where available 
does not alter the stellar masses significantly and our mass 
functions are consistent i n the overlapping mass range with 
the ones derived by Ilbert et al.1 (l2009b) inclu ding th e IRACl 
channel. In a similar analysis Fontana et al.l (|200^ find that 
including IRAC 1 -IRAC4 at z < 1 data does reduce the uncer- 
tainties in the mass estimate somewhat, but does not change 
the masses systematically (Fig. 1 in their paper). The reason 
for this is that in the redshift range of interest here (z < 1), the 
restframe near-IR light is sampled well enough by the J and 
Ks data, so that the addition of IRACl data does not add much 
information about the presence of an old stellar population. 

In addition to the above magnitude limits, we require that 
a photometric redshift be assigned to the objects and that the 
object not be in a masked image region (e.g. near bright stars, 
image border, detector defects) in any of the B j, Vj, t, and z^ 
bands. The total area of the survey we use after masking is 
1.73 square degrees. 

We use two methods to remove stars from the sample. First, 
sources where the best-fitting stellar SED template has a lower 
than the best-fitting galaxy SED template are discarded. 
This method is discussed and shown to be very reliable at 



TABLE 1 
Sample Size 



Zmin 


Zmax 


Volume 


Number of Galaxies 






10* Mpc' 


Active 


Passive Total 


0.2 


0.4 


0.56 


27931 


3553 31484 


0.4 


0.6 


1.23 


27258 


2289 29547 


0.6 


0.8 


1.92 


34298 


2587 36885 


0.8 


1.0 


2.53 


37131 


2954 40085 



t < 24 mag (lllbert et al.ll2009at 2% of the total population 
at < 24 could be stars not recognized by the SED classifier 
while only 0.2% of extended sources in ACS are classified as 
stars by the x^-criterion). In addition, we remove sources with 
point-like ACS r norph ologies using the catalog presented in 
iLeauthaud et al.l ( 12007 ). This ACS cut removes an additional 
1-3% of the total number of sources and is especially im- 
portant among passive galaxies at z > 0.6. We also note that 

15% of point-like sou rces in the COSMO S field are found 
to be galaxies, however, iRobin etal] (l2007h show that these 
are faint high-z objects that will not affect our analysis. 

Our final sample contains 138001 sources in the redshift 
range 0.2 < z < 1. We will divide this redshift range into four 
equally spaced bins in the following analysis. Table[T]lists the 
number of sources in our final sample and its subsamples as 
well as the volume probed by each redshift bin. 

2. 1 . Reliability of the Photometric Redshifts 

While our mass functions yield consistently steep faint-end 
slopes in all of the redshift bins studied, the lowest redshift bin 
provides the deepest probe of the low-mass galaxy population 
and our strongest constraints at this mass range. A top con- 
cern, therefore, is that distant galaxies with much larger red- 
shifts may have been mistakenly shifted to the 0.2 < z < 0.4 
bin and, with faint apparent magnitudes, assigned masses with 
logM/M0 < 9. Indeed, photo-z estimates often suffer from 
degeneracies between the 4000A break at z 0.2 and the Ly- 
man break at z ~ 4, especially if f/-band or near-IR photome- 
try is shallow or missing. We largely mitigate these problems 
by requiring detections for all sources in i"*" and Ks in addition 
to implementing our magnitude limits to reject noise and false 
detections. We confirm this by verifying that our results do 
not change if we reject all sources which show a second peak 
in their redshift probabiHty distribution with ampUtude larger 
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than 5%. We note that lllbert et alj (l2009al ) find no evidence of 
redshift degeneracies or persistent catastrophic outliers in the 
COSMOS photo-z catalog at our magnitude limits that could 
threaten the robustness of our result. 

Beyond the photometry requirements, we have searched the 
NUV-R-J color-color diagram for sources with unusual colors 
that could indicate a contaminating population. This is illus- 
trated in Fig. [T] Given the mass completeness limit estimate 
for each redshift bin, Mu^ ,, we highlight in blue those galax- 
ies with M < 3 X Miiin. These represent the lowest masses at 
which the star-forming sample is still complete. In large part, 
they follow the locus of star-forming galaxies. No outlier pop- 
ulation is apparent. We have also repeated this experiment for 
galaxies whose photo-z uncertainties are particularly large (a 
68% confidence interval larger than 0.5 in redshift) and again 
find no evidence of contamination. Finally, we use Monte- 
Carlo simulations described in more detail below to convince 
ourselves that the mass function is sufficiently robust against 
the errors in the photometric redshifts. 

2.2. Separating Star-forming from Passive Populations 

The low-mass galaxies contributing to the steep faint-end 
slope of the mass function are too faint to be morphologically 
classified using the available HST data. Instead, we turn to 
their colors and SEDs to shed light on their properties. Mo- 
tivated by the ability to break degeneracies between age and 
dust, we use a UV-optical-IR (NUV-R-J) color-color diagram 
to study our sample, finding that the best-fitting SED tem- 
plates returned by the photo-z code do a good job of distin- 
guishing star-forming from passive galaxies. We group to- 
gether SED types 1-8 fr om the COSMOS phot o-z catalog as 
"SEP early-types" (see l llbert et al.ll2009al and iPoUetta et al.1 
120071 for details on the SED templates) and highlight them in 
red in Fig.[T] At most redshifts, the early-type SED designa- 
tion overlaps well with the associated clump of passive galax- 
ies identified with re d N UV-R colors and blue R-J colors (see 
IWilliams et al.ll2009l and lllbert et al.ll2009 a). A joint color-cut 
indicated by the two lines would select additional passive sys- 
tems that are not selected by the SED classification. At z < 0.6 
the color cuts select ^ 10% more passive galaxies, many of 
these systems scattered between star-forming and passive se- 
quences. At z 1, however, the SED method identifies as 
many as 15% fewer galaxies in the passive clump. 

With none of these classifiers being demonstrably prefer- 
able, we choose to proceed with the SED-based classification. 
We have re-analyzed the data with the NUV-R-J-based classi- 
fication, and find that the results are fully consistent. In partic- 
ular, this excludes that the mass function of faint red galaxies 
is distorted significantly by some faint blue galaxies (which 
are more numerous by a factor > 10) being misidentified as 
red galaxies by one of these classification methods but not the 
other Also, if misclassified blue galaxies were (partly) re- 
sponsible for the faint red population, we would expect this 
population to become more numerous at high z where S/N is 
lower This turns out not to be the case (see below). 

In Fig. 12] we show the distribution of the galaxies in our 
sample in the absolute magnitude, Mr, versus u-r restframe 
color plane. The greyscale shading marks the density of all 
objects in this plane, while the blue and red contours outline 
the distributions of the blue and red subsamples, respectively. 

3. DERIVING STELLAR MASSES 

Stellar masses (the mass locked up in stars) are routinely de- 
rived by comparing multi-band photometry to a grid of stellar 



population models of varying star formation histories (SFHs), 
ages, and dust content. We follow the same approach in this 
work. Our procedure is described in detail and compared 
against spectroscopic and dynamical mass estimates in Drory, 
Bender, & Hopp ( 2004), although in the present work we will 
explore a wider set of stellar population models described in 
the following paragraphs. 

3.1. Stellar Population Models 

We allow for SFHs consisting of two components: a main 
population of stars formed in a smooth SFH and a modula- 
tion by a second population formed in a burst of star forma- 
tion just prior to observation. The main component is pa- 
rameterized by a SFH of the form ip(t) oc exp(-f/T), with 
r e [0.5, oo] Gyr and a metallicity of -0.6 < [Fe/H] < 0.3. 
The age, f, (defined as the time since the onset of star for- 
mation) is allowed to vary between 0.5 Gyr and the age of 
the universe at the object's redshift. We linearly combine this 
with a starburst modeled as a 100 Myr-old constant star for- 
mation rate episode of solar metallicity. We restrict the burst 
fraction, (3, to the range < /3 < 0.15 in mass. Higher val- 
ues of (3 are degenerate and unnecessary since such SEDs 
are covered by models with a young main component. We 
adopt a Chabrier (2003) IMF truncated at 0.1 and 100 Mq 
for both components. This choice of IMF is described by a 
power-law at M > 1 Mq and a log-normal distribution be- 
low. The Chabrier IMF yields masses lower by about 0.25 dex 
compared to a Salpeter (119551) single-slope IMF of the form 
dn/dm oc m~^'^^. We use the stellar population synthesis 
codes of lBruzual feCh arlot (2003, BC03) and the BC03 mod- 
els updated in 2007' with an improved treatment of TP-AGB 
stars, hereafter referred to as BC07 models. 

Additionally, both components are allowed to exhibit a dif- 
f erent and varia ble amount of extinction by dust. We assume 
a lCalzetti et alj ([2000) form for the extinction law. This takes 
into account the fact that young stars are found in dusty envi- 
ronments and that the starlight from the galaxy as a whole may 
be r eddened by a (g eometry dependent) different amount. In 
fact. lStasiriska et aTl (f2004) find that the B aimer decrement in 
the SDSS sample is independent of inclination, w hich, on the 
other hand, is driving global extinction (see, e.g..' Tullv et alj 
119981) . We verified that using the Milky Way extinction law or 
the SMC extinction law instead does not change our conclu- 
sions. This is due to our restricted redshift range which means 
we do not probe too far into the restframe UV where the dif- 
ferent slopes of the extinction laws matter most and because 
low-mass galaxies at z < 1 are not heavily obscured, unlike 
galaxies of the same mass at z > 2. 

We compute the fuU likelihood distribution on a grid in 
this 6-dimensional parameter space (t, [Fe/H], t, Ay ,(3, Ay), 
the likelihood of each model being cx exp(-x^/2). To com- 
pute the likelihood distribution of M/L, we weight the M/L 
of each model by its likelihood and marginalize over the stel- 
lar population parameters. The uncertainty in M/L is obtained 
from the width of this distribution, which falls between ±0.1 
and ±0.3 dex at 68% confidence level on average. The un- 
certainty has a weak dependence on mass (increasing with 
lower S/N photometry), and much of the variation is in spec- 
tral type: early-type galaxies have more tightly constrained 
masses than late types. 

We have also compared the masses obtained with this code 

' Often referenced as lCharlot & Bruzuall <2007h 
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Mr 

Fig. 2. — The distribution of galaxies in our sample in tiie absolute magnitude, Mr, versus u—r restframe color plane. The greyscale shading marks the density 
of all objects in this plane, while the blue and red contours outline the distributions of the blue (active) and red (passive) subsamples, respectively. 



to the ones obt ained with the code employed by Bun dy et alJ 
(I2006h and the iBruzual & Charlod (.2003) stellar population 
synthesis models. We find that the masses computed by both 
codes are consistent and show no trends with luminosity or 
redshift. Hence the stellar masses are not sensitive to the par- 
ticular choice of grid parameters and fitting procedure (within 
reasonable limits, of course). 

Of much graver concern than errors stemming from uncer- 
tainties in the photometry are systematic problems with the 
stellar population libraries and the implicit priors on evolu- 
tionary histories introduced by the finite model grid. For a 
comprehensive analysis see Marchesini et al. (2008). One of 
the main issues is the treatment of the thermally pulsing hor- 
izontal giant branch (TP-AGB) stars, whose light dominates 
in intermediate-age ( 0.5-2 Gyr) populations from the /-ban d 
through the near-IR (lMarastonll200l iMaraston et"aLl 120061) . 
The contribution of these stars can lead to changes of up to a 
factor of two to three in the stellar mass estimates for extreme 
cases (however, these are unlikely to be observed at < z < 1) 
and if the relevant wavelength range is poorly covered by the 
data. TP-AGB spectral signatures are strongest at ~ 1 /im 
rest-frame, which is always within the observed-frame colors 
in our dataset. Nevertheless, we base our model grid on the 
BC03 models as well as on the BC07 models which incorpo- 
rate an improve d prescription for the treatment of the TP-AGB 
phase (see also IBruzualll2007h . 

Despite the probl ems with late stellar evolutionary phases, 
IConrovet al.l(l2009l) in a thorough investigation of population 
synthesis modeling, argue that the M/L ratios are largely re- 
sistant to the uncertain contribution from TP-AGB stars as 
well as other limitations in the models. We will nevertheless 
compare stellar masses based on both evolutionary synthesis 
codes to make sure our results do not depend on a specific 
choice of a stellar population model. 

3.2. Mass Completeness Limits 

Determining the selection function in stellar mass for a 
flux-limited sample is not possible in a rigorous way: knowl- 
edge of the intrinsic frequency distribution of mass-to-light 
ratios of the galaxy population above and below the flux 
limit is necessary, but generally not available. Several meth- 
ods to (conservatively) estimate the completeness limit have 
been employed. Most commonly and simplistically, one as- 
sumes that a maximally old population (and hence with max- 
imally high M/L in the absence of dust) at the flux limit 
of the survey will yield a conservative upper limit to the 



TABLE 2 
Mass completeness limits. 







Active 


Passive 


^min 


^max 


logMu,„/M0 


logM,in,/M0 


0.2 


0.4 


8.3 


8.9 


0.4 


0.6 


8.9 


9.2 


0.6 


0.8 


9.2 


9.8 


0.8 


1.0 


9.4 


10.1 



mass that potentially could be affected by incompleteness 
(e.g. Drorv et al. 2003, 2005; Fontana et al. 2006; Borc h et alj 
|2006; Bundv et al. 2006; P erez-Gonzalez et aL.2008.) . 

Marchesini et al. (2008) employed a different approach: 
they use successively deeper datasets to empirically deter- 
mine their mass completeness function at shallower flux lim- 
its. This approach, while clearly better than the maximum- 
M /L method is not feasible for large, uniformly deep surveys 
such as COSMOS (there are no significantly deeper datasets 
observed in the same passbands). However, it is possible in 
a limited way using smaller pencil-beam surveys and we will 
do so as a sanity check. 

We wish to separate our sample into blue, star-forming 
galaxies and red, quiescent ones. It is therefore also nec- 
essary to determine separate completeness limits for both 
populations. We will refrain from attempting to correct 
for incompleteness. Instead, we conservatively estimate 
the point from which incompleteness may significantly af- 
fect our sample and limit ourselves to masses larger than 
this number. For the red population, a maximal M/L ap- 
proach seems sufficient. At the redshifts of interest (0 < 
z < 1), this population consists of ellipticals and quies- 
cent spirals with little evolution in number density with 
masses logM/Mp > 10 (Bell et al. 2004; Bundv et al. 20051 
Pozzetti et al. 2 007; .Perez - Gonzalez et al. 2008; Ilbert et al] 
2009b; Willia ms^al.1 120091) . Stronger evolution and dust- 
extincted starbursts appear in large num bers at z > 1.5 
(Fontana et al. 2006; Daddi et al. 2007; Perez -Gonzalez et"an 
l2008t illbert et al.i i200 9b). In the local universe most low- 
mass passive galaxies are dwarf spheroidals which are very 
unlikely to contain significant amounts of dust. We obtain an 
estimate for the completeness limit of our sample at < 25 
and Ks < 24 of logM/Mg = 9 at z = 0.3 and logM/M© = 10.1 
atz = 0.9. 

Blue star forming galaxies have much more varied M /L val- 
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ues. However, the largest variation is in ^ L* spirals, with a 
large contribution to this variation coming from variable dust 
extinction; blue objects fainter than M,- ^ -20 at z ~ have 
much more uniform properties (Kauffr nann et alJ 120031) and 
in generally much lower amounts of dust (Ahq < 1 mag at 
log M/Min < 9 com pared to Ana ~ 2 ± 1.5 at logM/M© > 
10; iBrinchmann et a l. 2004). In accordan ce with the con- 
clusions reached bv^ Kauffmann et al.l (120031) in the local uni- 
verse, the mass-to-light ratios at magnitudes fainter than Mr ^ 
-20 become roughly independent of magnitude, converging 
to M/Lr ^ 1. At brighter magnitudes than Mr ^ -20, M/L,. 
rapidly increases with luminosity. 

We can therefore assume that models with M/L^ ^ 1 resem- 
ble the faint galaxy population and a limiting mass computed 
using a range of models with values of M/L, close to or above 
1 will yield a sufficiently robust limit. Using plausible num- 
bers for a rather extreme model of a low-mass galaxy with 
a (light-weighted) age of 1 - 2 Gyr, star formation timescale 
(r = 5 - 10 Gyr), and the maximum dust extinction found in 
the local universe for such objects {Ana = 1 mag), therefore 
pushing M/L to the high side of the distribution, we obtain 
completeness limits roughly 0.4 to 0.7 dex lower in mass than 
for the red galaxies. We can verify these limits using th e much 
small er but deeper FORS Deep Field (FDF) data ( Drory et al l 
120051) . Comparing the COSMOS mass function at 0.4 < z < 
0.6 with the FDF mass function bin 0.25 < z < 0.75 from 
iDrorv et all (1200 5'), we see that the point where the COSMOS 
data start to fall below the FDF data coincides with the mass 
limit for the blue population as estimated here (see Fig.|5]be- 
low). The completeness limits are summarized in Table|2] 

4. THE GALAXY STELLAR MASS FUNCTION 

We compute the galaxy stellar mass function using the Vm^x 
method to account for the fact that fainter galaxies may not 
be detectable throughout a whole redshift bin. Hence, each 
galaxy in a given redshift bin contributes to the number den- 
sity an amount inversely proportional to the fractional volume 
of the bin in which the galaxy would be visible. Given the 
conservative mass completeness limits applied in our analy- 
sis, we find that the Vmax -corrections to the mass function data 
points above the mass limit are negligible (below 5% in all 
bins). Therefore, the data used in the analysis of the mass 
function are essentially free of incompleteness corrections. 

4. 1 . Uncertainties due to Photometric Redshifts 

Aside from the potential for catastrophic failures (discussed 
in § 12. It , the use of photometric redshifts implies imprecise 
distances and therefore imprecise absolute magnitudes or stel- 
lar masses. On the exponential part of the mass function, un- 
certainties in the distance will scatter objects preferentially 
from lower masses (where objects are exponentially more 
abundant) to higher masses. It may also introduce a system- 
atic effect on the faint-end slope. 

We investigate the effect of uncertainties in the photometric 
redshifts on the mass function by means of Monte-Carlo sim- 
ulations. We re-assign redshifts to each object drawing from 
the redshift probability distribution in the 95% confidence re- 
gion around the main peak. We take the asymmetry of the 
redshift probability distribution into account by using differ- 
ent confidence limits for the region below and above the main 
peak. We then recompute stellar masses with such new red- 
shifts and redetermine the mass function. In Fig. |3] we plot 
50 realizations of mass functions obtained from these sim- 
ulations. We conclude that redshift errors do not add large 




log M* [hyo"^ Msun] 

Fig. 3. — Monte-Carlo realizations of the stellar mass function to explore 
the influence of the uncertainties in the photometric redshifts on the shape of 
the stellar mass function. The mass function of passive galaxies is marked 
in red, that of star-forming galaxies in blue, and that of all galaxies in black. 
Each line con'esponds to one realization where the redshifts have been ran- 
domly drawn for each object individually from the 95% confidence region of 
the redshift probability distribution. 




log M* [hyo"^ Msun] 

Fig. 4. — The z ~ galaxy stellar mass function from SDSS data 
(Baldrv et alj|2008l ; solid line) convolved with the mass error distribution 
due to the use of photometric redshift based distances (dashed lines) at 
zG {0.3,0.5,0.7,0.9}. 



uncertainties at the faint end. The uncertainties in each mass 
bin at logM/M© < 10.5 we derive from these simulations 
are < 0.03 dex for the blue and total mass functions, and 
< 0.05 dex for the red mass function. The inferred faint-end 
slopes agree to within a fraction of the formal fit uncertainty 
on the slope parameter. 

The slope may still change systematically, though. To test 
this, we run another set of simulations, this time using only 
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objects with small redshift errors (the best 20% of the Az/(l + 
z) distribution) and increase their errors by a factor of 2 and 
3. In neither case is the slope significantly altered compared 
to the one obtained with the original small redshift errors (the 
objects in these simulations are obviously not quite as faint 
as those in the whole sample, but faint enough below M* to 
securely measure a change in slope). 

However, as expected, there is an effect on the bright end of 
th e mass function. Fi g.Hshows the z^O SDSS mass function 
by lBaldry et al.l (|2008) and the result of convolving it with our 
mass error distributions. The presence of mass uncertainties 
leads to a significant overestimate of the mass function at the 
high-mass end. In summary, we conclude that redshift uncer- 
tainties do not influence the faint-end slope of the mass func- 
tion systematically but do add a small random uncertainty. 
At the bright end they lead to larger uncertainties and a sys- 
tematic "brightening" of the exponential cutoff. We will add 
the uncertainties from these simulations (the scatter seen in 
Fig.O to the total uncertainties in our mass functions. 

4.2. Characterizing the Shape of the Mass Function 

We plot the galaxy stellar mass function in Fig.|5] The mass 
function of passive galaxies is shown in red, that of active 
galaxies in blue, and that of all galaxies in black. Filled sym- 
bols mark data points above the formal completeness limit, 
open symbols data below the completeness limit and hence 
where number densities are likely to be artificially low. 

It is evident from Fig. |5] that a single Schechter function 
would not be a good fit to the total mass function (black 
points) or to the active galaxy mass function (blue points) at 
any redshift z < 1. At z < 0.5, the passive galaxy mass func- 
tion clearly shows a marked paucity of objects at intermediate 
mass and a sharp upturn at low mass. At z 0.7 the upturn in 
the passive galaxy mass function is still detected, however the 
mass function is no longer complete at any point that deviates 
upwards, and at z ~ 0.9 we are unable to detect an upturn. 
At these two redshift bins, using only data above the com- 
pleteness limit, the passive-galaxy data are adequately fit by a 
single Schechter function. 

One way to describe the mass function shape is as a sum 
of (at least) two components, that is, a bimodal distribution. 
Double Schechter functions (the sum of two Schechter func- 
tions, sometimes sharing a common normalization, </)*, char- 
acteristic scale, M*, or both) have been used by other au- 
thors to fit the luminosity or mass function s with steepen- 
ing faint-end part s (jPriver et al. 1994; Pop esso et alJ l2006l 
iBaldry et alJl2008t IPozzetti et al.,,2009,) . These work well to 
describe the total population and the active sub-population in 
our data set. However, more flexibility is needed to model 
the transition region between bright and faint populations in 
the passive sub-population at z < 0.6, as these functions fail 
to reproduce the fairly wide dip in the data. We therefore opt 
to fit all mass functions by a sum of two Schechter functions 
without restricting the parameters further: 



(j>(M)dM= (j)i,(M)dM+ (j)f(M)dM ■■ 



M 

M* 
M 
M* 



exp 



exp 



M , 
dM+ 

M , 
dM. 



(1) 



We define and such that > thereby identifying 
the first term, (l)b(M), with a population of bright galaxies and 



the second term, <f>f(M), with a population of faint galaxies. 

The results of fitting Eq.[T]to the data are shown in Fig.|5] 
Solid lines mark the fit, 0(A/), to the active, passive, and to- 
tal mass functions while the dashed lines show the individual 
components, (t>h{M) and <j)f{M), of each fit. We again use blue 
for the active galaxies, red for the passive galaxies, and black 
for the total population. The best-fit parameters of Eq. [T] are 
listed in Table [3] The table also lists the reduced values 
for the best fit, as well as the obtained by fitting only a 
single Schechter function to the data instead of Eq. [T] The 
reduced x^ values are between 1 and 2. The best x^ values 
obtained with a single Schechter function are between 7 and 
13, confirming that a single Schechter function indeed does 
not provide satisfactory fits with high significance. 

We find a faint-end slope significantly steeper than previous 
studies, af ^ -1.7 ±0.15 with very little variation with red- 
shift and between the blue, red, and total samples. With shal- 
lower mass limits, slopes typically a round a ~ -1 .2 are found , 
(e.g . Bell & de Jon2 2001; Cole et al'l200ll:lDrory et alJ2004l: 
iFontana et alJ 12004 i2006c illbert et al.ll2009U although evi- 
dence for a steeper slope has been put forward, for exam- 
ple bv Drorv & Alvarea (120081) from comparing the evolution 
of the mass function to the expectation from star formation 
rates. The value we find at z=0.3, a t ^ -1.7±0.1, is consis- 
tent with that of Baldry et alJ (12008b in the SDSS, who find a 
best-fit value of a2 = -1 .58 ± 0.02, with 0:2 = -1 -8 providing 
equally good fits given the non-negligible role of systematic 
surface brightness selection effects at their faint end. We find 
that the faint-end slope of the blue and the red population at 
z = 0.3 and z = 0.5, where both can be measured, are remark- 
ably similar. The red faint-end slope is formally steeper by 
0.1, but with a large uncertainty of 0.2-0.3. (Evidence for 
an excess of red dwar f galaxies was recently also found by 
'Salim beni et"aLl (120081) in the GOODS data). We do not detect 
a significant trend of the faint-end slope of the blue population 
and hence of the total mass function with redshift up to z = 1 . 

The characteristic mass of the faint blue sub-component 
shows a slight increase with redshift, from logMj 9.6 at 
z = 0.3 to logMf* ^ 9.8 at z = 0.9 (significant at the 2-a level), 
while the characteristic mass in the bright component, 
does not show evolution in this dataset. 

Surveys of 1 - 2 deg^ based on single fields can suffer 
from significant cosmic variance uncertainties. The volume 
probed by our redshift bins lies between 0.56 x 10* /i^g Mpc-* 
and 2.53 x 10* /lyg Mpc-' (equivalent to box sizes between 
82 hjQ Mpc and 136 /i^q Mpc; see Table [T]). Cosmic vari- 
ance on such scales is mostly thought to affect the overall nor- 
malizati on of the ma ss function, although possibly the shape 
as well dStringer et alJl2009l) . In the COSMOS data, we no- 
tice a str ong underdensity at z ~ 0.5, as has been seen before 
(see also lMeneux et al.ll2009l) . with the characteristic density, 
(p*, being about a factor of two lower (more so in the red 
population than in the blue). Looking at the dist ribution of 
spectr oscopic redshifts from zCOSMOS (Fig. 8 in lLillv et alj 
I2009h . it is evident that over the interval 0.4 < z < 0.65 there 
is a marked paucity of dense structures compared to lower 
and higher redshifts, which confirms our above observation. 
The comparison with the spectroscopic zCOSMOS sample 
also excludes the possibility that the density structure in the 
field is caused by systematic failures of the photometric red- 
shifts. This variance in average density influences our ability 
to make statements on redshift evolution significantly. De- 
spite the large cosmic variance, we do see a shift from blue to 
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log M* [h7o~^ Msun] 

Fig. 5. — The stellar mass function of galaxies in four redshift bins in the interval 0.2 < z < 1.0. The mass function of passive galaxies is marked in red, that 
of stai-forming galaxies in blue, and that of all galaxies in black. Data points below the mass completeness limits ai'e denoted by colored open symbols. The 
uncertainty in the mass function due to Poissonian err ors in the counts as well as the uncertainty propagated through photometric redshifts (see text) is shown 
as shaded regions. Data from the FORS Deep Field i Drorv etalj|200^ are shown as small black open squares at z = 0.5 for comparison. The solid lines show 
double Schechter function fits to the data (Eq.[T). The thin dashed lines show the individual bright and faint components of the double Schechter functions. 



red galaxies at the bright end as redshift decreases. 

A word of caution is warranted at this point; Eq.[T]is by no 
means the only possible fitting function able to reproduce the 
data in a statistically satisfactory way. For instance, a faint 
component with a sharper than exponential cut-off works just 
as well. Also, the bright part of the (red galaxy) mass func- 
tion can be fit by a symmetrically peaked function, similar 
to a Gaussian. Thus, one should not base any interpretations 
solely on the particular fitting function chosen. We empha- 
size that the fitting functions in the overlap region between the 
bright and faint mass functions should not be taken to be more 
than descriptions of the data. Without classifying objects in- 
dividually as belonging to the bright or faint subpopulation, 
the transition region between the two can be fit in many ways 
and the particular representation chosen here should not be 
assigned physical meaning. 

4.3. Effects of Stellar Population Models 

Can the features of interest in the stellar mass function be 
explained as artifacts introduced by the adopted stellar popu- 



lation models? First we test whether TP-AGB stars can influ- 
ence the faint-end slope of the blue star-forming mass func- 
tion. We recompute the mass function using a stellar popula- 
tion grid based on the BC07 mo dels, keeping all other model 
grid parameters identical (see § 13.11 ). Fig. |6] shows the dif- 
ference in the mass function obtained with the BC03 models 
(Fig- 13 and the mass function obtained with the BC07 mod- 
els. We do not find significant systematic differences in the 
mass function of passive galaxies. In contrast, systematic dif- 
ferences are small but noticeable in the active population and 
they increase with redshift. As expected, the stellar masses 
obtained with the BC07 models are smaller than the on es ob- 
tained with flie BC03 models (^araston 2005; Marasto n~et alJ 
'2006';' Bruzuaill2007h . The difference is negUgible at z - 0.3. 
Beyond z ~ 0.5, the difference exceeds the uncertainties in 
the masses over much of the mass range, increasing at masses 
above and below lO'^M© to a maximum of 0.1 to 0.15 dex. 
However, these differences are not large enough to change the 
shape of the mass function and the conclusions reached in this 
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TABLE 3 

DOUBLE-SCHECHTERFIT PARAMETERS TO THE MASS FUNCTION 



z 






10 hjQ Mpc dex 


logM* 




10 hj^Mpc dex 


logMj 


Of 


, 2 
X 


2 a 
^single 


0.2. . 


.0.4 


All 


2.89 ±0.23 


10.90±0.11 


-1.06 ±0.03 


1.80 ±0.29 


9.63 ±0.09 


-1.73 ±0.09 


1.9 


13.0 






Passive 


1.96 ±0.12 


10.80 ±0.99 


-0.49 ±0.14 


0.49 ±0.1 2 


9.54 ±0.09 


-1.85 ±0.20 


1.8 


16.8 






Active 


1.66±0.10 


10.83±0.10 


-1.23 ±0.03 


0.98 ±0.23 


9.66±0.10 


-1.75±0.12 


1.8 


7.0 


0.4. . 


.0.6 


All 


1.74 ±0.09 


10.91±0.11 


-1.05 ±0.02 


1.43 ±0.23 


9.70±0.10 


-1.76±0.16 


1.6 


13.6 






Passive 


0.95 ±0.09 


10.85±0.11 


-0.40 ± 0.06 


0.28 ±0.03 


9.41 ±0.19 


-1.84 ±0.29 


1.6 


12.7 






Active 


1.38 ±0.06 


10.76±0.11 


-1.13±0.05 


1.29 ±0.04 


9.69 ±0.09 


-1.71±0.19 


1.3 


9.6 


0.6. . 


.0.8 


All 


2.16±0.13 


10.95±0.10 


-0.93 ±0.04 


2.89 ±0.26 


9.75±0.10 


-1.65 ±0.08 


1.5 


11.9 






Passive 


0.90 ±0.05 


10.94 ±0.09 


-0.39 ±0.05 










0.7 






Active 


1.86±0.11 


10.80±0.10 


-0.95 ±0.04 


2.78 ±0.31 


9.75±0.10 


-1.61±0.11 


1.6 


9.1 


0.8. . 


.1.0 


All 


2.94 ±0.1 3 


10.92±0.10 


-0.91 ±0.03 


2.12±0.29 


9.85±0.10 


-1.65 ±0.24 


1.4 


9.1 






Passive 


1.03 ±0.05 


10.91±0.11 


-0.29 ± 0.04 










0.3 






Active 


2.51±0.10 


10.81±0.11 


-0.97 ± 0.03 


2.16±0.37 


9.80±0.10 


-1.66 ±0.36 


1.4 


8.0 



Reduced x~ for ^ fit with a single Schechter function 
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log M* 

Fig. 6. — The ratio of the galaxy stellar mass function computed with stellai' population synthesis models based on the BC03 code and the BC07 code which 
includes a refined treatment of the post-AGB stellar evolutionary phases. The values for active and passive galaxies are marked by blue and red symbols, 
respectively. The shaded regions denote the statistical uncertainty in the mass functions. 
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Fig. 7. — The luminosity function of galaxies (the number density of galaxies per comoving volume) in the restframe B, r, i, and z bands at z = 0.5. Passive 
galaxies are marked in red, star-forming galaxies in blue, and all galaxies in black. Data points below the completeness limit are denoted by open symbols. 



work in any significan t way (see also IConroy et al.ll2009l and 
lMarchesinieral]l2008l) . 

Could the mass function shape be affected by a failure to 
correctly model the very low M/L population? We plot the 
luminosity function of galaxies in the rest-frame B, r, i, and 
z bands at z ~ 0.5 in Fig. |7] The B-band LF shows no bump 
in the blue or total galaxy populations; however a bump-like 
feature appears in the LF of blue and of all galaxies as one 



moves to redder filters, and is unmistakable in the z-band LF. 
So the dip and steep faint-end slope are present in luminosity 
functions redwards of the restframe B band and do not arise 
from converting luminosity to stellar mass via SED fitting. 

5. DISCUSSION 
5.1. Comparison to Previous Results 
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log M* [hyo"^ Msun] 

Fig. 8. — A comparison of the mass function deduced in this work with data from the hterature. The mass function of passive galaxies is marked in red, that 
of star-forming galaxies in blue, and that of all galaxies in black. The solid lines show double Schechter function fits to the data (Eq. [T] the thin hues show 
the individual bright and faint components of the double Schechter function). The fit at low z is repeated in higher-^ panels as dotted Hues. The magenta line 
denotes the Baldry et al. (2008) z ~ SDSS-based mass function convolved with the photo-z error distributions (see text). The dashed blue, red, and black lines 
are the 3.6/^m-selected mass functions of active, passive, and all galaxies, respectively, in the COSMOS survey by Ilbert et al. (2009b). The green dash-dotted 
line is the mass func tion by Perez-Gonzalez et aU >2008i) . Additionally, we plot the Press-Schechter halo mass function scaled by the global baryon fraction /j, 
IDunklev et aU2009l ; WMAP5) in cyan. 



Before discussing possible interpretations of the more com- 
plicated mass function shape apparent in the COSMOS data, 
it is useful to compare our results to previous work with 
the aim of determining whether evidence for similar behav- 
ior has been found in other surveys. In Fig. [8] we compare 
our mass function to a number of results from the literature 
and for later refe rence we include the halo mass function 
(iReed et alJ l2007h where the halo masses have been multi- 
plied by the global bar yon fraction, fh = ilb/ ^m, taken from 
the WMAP 5-yr data dDunMey et alJl2009). We show the 
z SDSS-based mass function by (Baldrv et al. 2008) con- 
volved with the ph oto-z error distributions as shown in Fig.|4] 
and discussed in § 14.11 We also plot the mass functions in 
the redshift range < z < 1 of active, passive, and all galax- 
ies in the COSMOS survey selected at 3.6/im by Ilbert et al. 
(t2009b) using the same photometric redshifts as the present 
w ork. Finally the total mass fu nction selected at 3.6-4.5/im 
bv lPerez-Gonzalez et al.l ( 120081) is also plotted. 



Where they overlap in mass, our mass functions agree very 
well with previous work, with some systematic differences 
stemming from differences in the underlying stellar popula- 
tion syntiiesisjibrariesjnd gri ds of star form ation histories 
(see i Marchesini et al.ll2008l and lConrov et al.li200 9 for a sys- 
tematic study on the influence of stellar population grid pa- 
rameters on mass determinations). Of notable importance, 
though, is the fact that the bright components in our two- 
component fit to the mass function agre e very w ell with the 
blue and the red sub-components in Ilb ert et al.l (|2009 b), de- 
spite (minor) differences in the definition of blue and red 
galaxies. This gives us confidence that the decomposition 
of the mass function into a bright and a faint component 
with six free parameters is not strongly degenerate. Specif- 
ically, both the mass scale of the bright components and 
their faint-end extrapolations are compatible with the sin- 
gle Schechter function fits to the more restricted data sets. 
We also confirm the build-up of the faint part of the red se- 
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quence at z < 1 observed by other groups (e.g. iB ell et all 
2004; Bundvetal. 2006; Bell et al. 2007; Faber et al T]2007l; 
Perez-Gonzalez et alj2008i;Jlbert et al.i2009bti Williams et alj 



2009h . 



As we discuss further below, the apparent multi-component 
na ture of the ma s s func tion has been discussed a t low redshifts 
bv lBaldrvetai] (120081) and lLi& WhTtg (120091) and observed 
at some level in other high-z spectroscopic and hence signifi- 
cantly shallower studies of th e COSMOS field dPozzetti et alj 
l2009t iBolzonella etal]|2009h . The advantage of the current 
analysis is that it probes to lower mass at redshifts beyond 
z = 0.2, thus allowing a larger dynamic range which for the 
first time at high-z reveals more complicated behavior, specif- 
ically the steepening of the faint-end slope a nd an additional 
popul ation of faint red galaxies (see also ISalimbeni et alj 
120081) . A similar steepening in the luminosity function be- 
low M, ~ -17 has been convincingly detected recently in clus- 
ters (Driver et al. 1994; Trentham & Tullv 2002; Hilk er et al. 
2003; Popesso et al. 2005, 2006), groups fc entham & Tullv 
2002; Trentham et al. 2005; Gonzalez et aLP2006), and in the 
field tBlanton et al. 2005). Baldry et al. (2008) find that the 
local galaxy stellar mass function steepens as well below 
logM/M© - 9.5 (see also Salucci & Persic 1999). 

5.2. Halo Mass to Stellar Mass Relation 



o 



o 




log Mh 

Fig. 9. — Stellar mass as a function of halo mass (bottom panel) and ratio 
between stellar mass and halo mass (top panel) at redshifts of z = 0.3 and 
z = 0.9 determined by abundance matching the galaxy stellar mass function 
to the halo mass function. The relation determined by .Moster et aU C2009,) 
are shown for comparison. 



For the following discussion of the physical significance of 
the various morphological features of the mass function, it is 
useful to first place the mass function in the context of the halo 
mass function. We will do so by analyzing the halo mass vs. 
galaxy stellar mass relation, noting that the more complicated 
shape of the stellar mass function also leads to a more com- 
plicated relation between galaxy stellar mass, Mg, and halo 
mass, Mh. 



The relationship between the halo mass function and the 
galaxy (stellar) mass function can be written as 



,(Mg,f) = 



dMh 



dM„ 



6h(Mh,f). 



(2) 



Such a one-to-one correspondence between halo mass and 
galaxy stellar mass at some time, Mh(Mg), can be found by 
requiring that the cumulative number density of halos above 
a given mass and galaxies above a corresponding stellar mass 
be equal (abundance matching). 



,(M)dM = 



b^{M)dM. 



(3) 



This approach has been shown to be sufficiently accurate to 
match the tw o-point correlation function of galaxies and halos 
dConrov et al. 2006; Moster et al. 2009). 

To obtain a more realistic picture of the halo masses 
of galaxies, we must include sub-halos in our abundance 
matching procedure, as halos (and hence their galaxies) may 
survive for considerable time. Analyzing hig h-resolution 
cosmological N-body simulations, iGiocoli et al.l (2008) and 
lAnpulo et al.l (12008) provide fitting formulae for the number 
of sub-halos per host halo per logarithmic in terval of sub-hal o 
to host mass ratio, dNjd\n(Ms/M) (see also lGao et alj2004t) . 
This quantity can be converted to a sub-halo mass function, 
0sub(A/s), which is the more convenient quantity for our pur- 
poses by changing variables to obtain dN^/dM^, multiplying 
by the abundance of host halos, and integrating over host halo 
mass from Ms to infinity: 



(Ms): 



dN, 
dM, 



dNi, 1 



dN, 



Ms 



t/Mh Ms dln{MjM) 



dM. (4) 



This sub-halo mass function can then be added to the distinct 
halo mass function to obtain the total mass function that we 
use to match to the abundance of galaxies: 



(/)(M) = (^distinct + 0sub. 



(5) 



We note that we use the mass of the sub-halo at infall 
time, not the mass of the sub-halo at the time of observa- 
tion for matching abundances. Most of the galaxy's proper- 
ties will have been set by the time infall occurs (and shaped 
by the potential of the halo at that time) and galaxy growth 
is likely to slow down or even stop as soon as the galaxy 
becomes a satellite. This approach is motivated by mod- 
els of galaxy formation that generically predict that stel- 
lar mass is tightly linked to the potential well in which 
the galaxy forms dKa uffmannjEalJliggl ICole et all [19941; 
ISomerville & Primacklil99 9). Sub-halos loose mass by tidal 
stripping, however their stars are more centrally concentrated 
and are not stripped until most of the dark matter has been 
lost. The relevant mass scale for sub-halos, therefore, is the 
virial mass at infall time. 

We proceed using a halo mass function consisting of the 
distinct halo mass fu nction byfe eed et al. (2007) and the sub- 
halo abundance by Giocoli et alJ i2008) processed through 
Eq.m In Fig.|9] we show the resulting relation between halo 
mass, Mh, and galaxy stellar mass, Mc,at redshifts z = 0.3 and 
z = 0.9 obtained by abundance matching to our total galaxy 
stellar mass function. For comp a rison, we also show the re- 
lation obtained by Moster et al.l ( 120091) by populating halos 
from N-body simulations with galaxies using semi-analytic 
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methods and the requirement that the stellar mass function be 
reproduced. 

Without accounting for substructure, Mg(Mh) evolves rela- 
tively uniformly at all masses. With substructure, much less 
change in the relation at low mass is seen; there is little evo- 
lution in Mg(Mh) from z = 0.3 to z = 0.9 at logMh/M© < 11 .5 
(corresponding to logMg/M© < 9.5, about the mass of the 
LMC). This is to be expected if the stellar content of low- 
mass halos is limited by feedback. There is evolution in 
the Mg(Mh)-relation at masses above M», where the stellar 
mass per halo mass increases by a factor of 1.23 from z = 0.9 
to z = 0.3. The peak in star formation efficiency occurs in 
logM/Mo - 12.2 halos and reaches 23% at z = 0.3 and 19% 
at z = 0.9; galaxy formation - if there is a monotonic relation 
between stellar mass and halo mass such as the one arrived at 
by abundance matching - is always an inefficient process. 

While the flattening of the Mg(Mh)-relation above M, (cor- 
responding to the exponential cut-off in the galaxy mass func- 
tion) is understood as an effect of the inefficiency of cool- 
ing in large halos (e.g.'R ees & Ostrike3ll977l; IWhite & ReesI 
119781) . the steepening of the Mg(Mh)-relation at halo mass of 
logM/MQ ~ 1 1 suggests another qualitative change in behav- 
ior of the feedback efficiency. 

The Mg(Mh)-relation is fundamental to the process of 
galaxy formation within dark matter halos and understanding 
what shapes this relation is equivalent to understanding the 
galaxy stellar mass function. At this point we wish to remark 
that, due to the derivative factor dMb/dM^ in Eq.|2] or equiv- 
alently, the integral nature of the abundance matching relation 
in Eq. |3] a single change of slope of Mg(Mh) leads to a dip- 
bump structure in the inferred galaxy mass function, given a 
power-law halo mass function. The more abruptly the change 
of slope occurs, and the larger it is, the more pronounced the 
dip-bump feature in the mass function becomes. The mass 
function shape we observe, therefore, leads to an approxi- 
mately double-power-law form of the Mg(Mh) relation, as can 
be seen in Fig.|9] 

5.3. Interpreting the Shape of the Mass Function 

There are several ways to interpret the shape of the total and 
type-dependent mass functions in COSMOS. Broadly speak- 
ing, we will concentrate on two interesting features, namely 
the "upturn" towards steeper faint-end slopes at low masses 
(logM/M0 ~ 9) and an apparent "dip" or plateau at masses 
just above logM/M0 ~ 10 that sets in below the traditional 
M* . We note that even breaking down the mass function into 
these two components is, in itself, an interpretation since the 
steepening faint-slope alone may provide a viable explana- 
tion. As we discuss below, our interest in separately identi- 
fying the dip is motivated by comparisons to the underlying 
halo mass function which is a strict power-law at the relevant 
masses. In this context, the dip represents a suppression of the 
stellar mass associated with halos in this range; this is a con- 
sequence of a change in dMh/dMg such that there is a smaller 
change in the mass of galaxies per change in the mass of ha- 
los. The purpose of this section is, hence, to explore several 
physical mechanisms that may explain both the upturn and 
the dip and to speculate on the role of future observations in 
distinguishing them. 

Usin g the zCOSMOS spectroscopic survey [Pozzetti et aH 
( 120091) interpret the mass function they recover as being com- 
posed of early-type galaxies dominating the massive part and 
late-type galaxies dominating the less massive part and con- 
tributing to the steep faint-end slope. They find that each of 



these components is well fi t by a Schechter functio n in the 
mass range they consider. iBolzonella et al.l (120091) use the 
same sample to investigate the bimodality as a function of 
environment. The find that at z < 0.5, the shape of the galaxy 
stellar mass function in high- and low-density environments 
become extremely different, with high density region s show- 
ing a stronger bimodality. Also, lllbert et ail (l2009bh remark 
that the 3.6/im-selected COSMOS mass function of all galax- 
ies is not well-fit by a single Schechter function. They prefer 
the sum of their (Schechter) fits to the red and blue galaxies to 
describe the whole population, as this reproduces the dip at in- 
termediate mass of logM ~ 10 and z < 0.6 they observe. Due 
to the IRAC selection, they do not reach faint enough limits 
to detect the steep faint-end slope or the bimodal structure at 
highe r z. We over plot the data from the PDF (iDrorv et al.l 
120051) at z = 0.5 in Fig.|5]a nd note that t he bim odality is visi- 
ble there as well, although IProry et aH (12005*) do not discuss 
it because a single Schechter fit seemed adequate at the time 
(apart from the data point at logM =11 which comes out too 
high). 

We emphasize that we see even more structure in the mass 
function than observed in the works discussed above; those 
studies only detect a dip in the total mass function and inter- 
pret its origin as being due to the superposition of active and 
passive populations which by themselves are unimodal. In 
this paper, we detect a significant change in slope in the mass 
functions of active galaxies as well as a marked bimodality in 
the mass function of passive galaxies. 

As mentioned above, previous work has attributed the two- 
component nature of the mass function to a superposition of a 
(unimodal) blue galaxy population modeled as a Schechter 
function and a red galaxy population, also modeled as a 
Schechter function but with a larger characteristic mass. In 
this picture, the bimodality arises through the transformation 
of blue galaxies into red galaxies in a process that must be 
linked to the mass evolution of a transitioning galaxy to de- 
plete the dip and create a bump at M ~ M* . In contrast, we 
find that the blue galaxy population by itself shows a dip sig- 
nature and a steepening of the power-law slope, which can 
be interpreted as arising from a bimodal distribution. In fact, 
both the red and the blue populations in Fig. |5]can be inter- 
preted as bimodal. However, the bimodality in the blue pop- 
ulation seems to be more pronounced at high z compared to 
low z: it becomes weaker as redshift decreases. This effect 
can also be seen in the values reported in Table [3] Indeed, 
the best values for single Schechter fits are obtained for the 
blue component at low redshift (Xsinoie ~ ^ O)- The signature 
of the bimodality (the dip at intermediate mass, or the bump 
at around M*) "moves" from the blue population at high z to 
the red population at low z. This is likely because some of 
the blue galaxies that make up the bump a round M* in the 
blue mass f unction t urn red with time (e.g. [ Bell et al.ll2004t 
iBundv et al. 2006; Fa ber et al.ll2()07HWilliams et al.ll2009l) . 

The fact that the blue mass function is itself bimodal at z ~ 1 
implies that the shape of the total mass function predates the 
emergence of the red sequence. We therefore suggest that 
the physical explanation must be a mass-dependent effect or 
mechanism that is largely separate from the process that trans- 
forms star-forming disks into passive spheroidals. In other 
words, this suggests a new dichotomy in the galaxy distri- 
bution besides the well-known red-blue distinction. We will 
follow this hypothesis in the following, using the behavior of 
blue galaxies to stand in for interpretations of the full popula- 
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tion and returning to red systems in the subsequent section. 

5.3.1. Blue Galaxies 

Before discussing possible physical explanations, we note 
that among star-forming, blue galaxies, there is an interest- 
ing separation between (giant) spirals of Hubble types Sa-Scd 
that dominate the mass function at the massive end and make 
up the majority of objects above the dip, and dwarf galaxies 
of Hubble types Sd-Sm/Im that dominate the less massive, 
power-law part of the blue galaxy mass function. Disk pa- 
rameters such as disk radius and luminosity remain roughly 
constant (with a large spread) along the Hubble sequence at 
types Sa-Scd. For types later than Sd galaxies, the Hubble se- 
quence turns into a luminosity sequence, with later type mor- 
phologies corresponding to sma ller and less massive gala xies 
(see, for example, the review by 'Robert s & Haynesll 1994 1 and 
references therein). It is also worth noting that Disk galax- 
ies without a significant bulge component exist (e.g., Hubble 
type Sd and later) and can be accurately characterized as pure 
disk systems, harboring at most a nuclear star cluster May we 
tentatively identify this distinction in disk properties with the 
two components in the blue population's mass function? If 
so, the physical processes discussed below must also account 
for these morphological differences. 

Star Formation Efficiency. The dip in the stellar mass 
function corresponds to a change in the stellar mass vs. halo 
mass relation such that there is a smaller change in the num- 
ber of galaxies when measured against the power-law shape 
of the halo mass abundance (Figs.[8]and|9]). It is possible that 
this deficit arises from a change in the efficiency of star forma- 
tion with mass. In this interpretation, the baryonic mass func- 
tion, if it could be observed, would show a close to power-law 
shape below M* , indicating a more smoothly varying fraction 
of Mbaiyon/A^haio with Mhalo- For this to still be the case at 
later cosmic times we must also require that the fraction of 
gas retained by any halo in the interesting mass range must 
be a smooth function of halo mass and also that this fraction 
not be too small. Then, a decrease in SFR efficiency below 
a given mass scale (say logM/M© ~ 10) would lead to an 
increasing gas fraction, /gas, at lower masses and a decreas- 
ing stellar mass fraction. This is an alogous to the break in 
the TuUy-Fisher relation dTully & Fisheilll977 ) where galax- 
ies with Vc < 90 km s"' fall below the relation defined by 
more massive galaxies. If plotted against total disk mass, 
Mgas +Mstar, instead of luminosity, a single relation is restored 
with Mgas +Ms,ar ^ ( McGaugh et al. 2000). The decrease in 
/gas as a function of M has been observed. Thus galaxies with 
baryonic masses near logM/M© ^ 10 would have observed 
M several factors lower, thereby creating a dip in the stellar 
mass function at logM/M© ^ 10, and leading to a steepening 
below logM/M© 9. 

Using a si mple fitting spl ine to observations of /gas(M) 
compiled by | Hppkins et a P (2009cj) including data from 
iBell & de JongI (12001b . lKannappanr (l2004h . and iMcGaughl 
(120051) . we transform the total stellar mass function at z = 0.3 
into an approximate baryonic mass function. The gas fraction 
hereby is defined as the ratio of gas mass to stellar plus gas 
mass, /gas = Mgas/ (Mgas +Mstai-). We plot the results in Fig.fTOl 

Using the observed gas fraction as a function of mass, 
the dip in the mass function can be lessened, but not re- 
moved. However, as expected, an appropriate assumption 
for fgas(M) leads to baryonic mass functions in which the 
dip can be "straightened" out (dot-dash line). However, the 
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Fig. 10. — The top panel shows the gas fraction, /gas, as a function of 
stellar mass. The data are taken from a compilation by Hopkins et al. 1 20093), 
including Bell & de Jong (2001, triangles) Kannappan (2004, circles), and 
IMcGaugh 12005, squares). A spline interpolation to the data is shown as 
a dashed line. The gas fraction necessary to remove the dip structure in the 
mass function is shown as a dash-dotted line. In the bottom panel, we plot the 
total stellar mass function atz = 0.3 and the baryonic mass function obtained 
by applying the gas fractions from the top panel. Power-law slopes of </> ~ 
M^- for the halo mass function and <j> ~ M"' ' for the stellar mass function 
are plotted for reference. 



amount of /gas(M) required is about a factor of 2 higher 
than current low-z observations (dashed line). This compar- 
ison involves several large uncertainties, but plausibly, the 
difference could be made up for in a warm phase that is 
not eas ily detected ( Cen & Osti-ike^ll999l; lOave et alj|2001l: 
ICen & Ostriker 2006). One would also expect /gas to increase 
with redshift in line with rising star formation rates. We are 
able to conclude that it is at least possible that the dip in the 
stellar mass function reflects a mass dependence in the ability 
for baryons to cool and form stars. 

What causes this mass-dependent SFR efficiency? We note 
that the deviation of the stellar mass function from power-law 
form at logM/M© 9- 10 is close to the scale where su- 
pernova feedback b e comes less irnporta nt, Vc ~ 100 km s"' 
(iDekel & Silk|[T986l: iBenson et al.ll2003l) . As long as super- 
nova feedback is efficient in regulating the conversion of gas 
into stars, the ratio of stellar mass to halo mass is set by the 
scaling relation of feedback efficiency with halo mass (or cir- 
cular velocity). Once the (DM) potential is deep enough that 
supernova feedback can no longer remove significant amounts 
of gas from the halo, the attainable ratio of stellar mass to dark 
matter mass becomes larger as a larger fraction of the gas in 
the halo can cool to form stars. As a consequence, the stel- 
lar mass function deviates upwards from its power-law form. 
(At the massive end, the stellar mass is limited by the cool- 
ing timescale becoming too large and the ratio of stellar mass 
to dark mass decreases again sharply causi ng the exponen- 
tial cut-off in the stel lar mass function (Rees & OstrikeJ 19771 : 
IWhite & Rees"1978). The transition from SFR efficiency sup- 
pressed by supernovae feedback to peak SFR efficiency may 
therefore lead to the observed dip in the stellar mass function. 

There are two potential tests for this scenario. The most ob- 
vious is to seek better measures of /gas as a function of M and 
redshift, but such observations remain challenging. A second 
test is to determine whether the mass scale that defines the 
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dip evolves with redshift. Because we expect /gas to be higher 
in more massive galaxies at early times (since at z > 1 many 
more are likely to be rapidly forming stars), the dip in the stel- 
lar mass function should move to higher masses at early times 
as well. Unfortunately, the mass functions presented here are 
likely to be too affected by cosmic variance to robustly mea- 
sure redshift evolution in the dip feature, but future data sets 
will be able to overcome this problem. Also, a more rigor- 
ous analysis taking into account cosmological gas accretion 
as well as a realistic feedback prescription and outflow model 
is needed, but is beyond the scope of this paper 

Hierarchical assembly. Another possibility is that the dip 
at logM ~ 10 and the bump around M* are formed by a de- 
pletion of stellar and gas mass at intermediate mass resulting 
from a pile-up of the end-products of mass-dependent assem- 
bly around M* . In this scenario, both the stellar and baryonic 
mass functions would exhibit dips at logM ~ 10 because the 
assembly of baryonic mass is accelerated at a certain scale. 

It is known tha t the halo-halo merge r rate depends very 
weakly on mass jFakhouri & Mai 120081) . However, if the 
stellar mass, or in this case, baryonic mass fraction, peaks 
in halos at a given mass scale — perhaps as a result of the 
SFR feedback discussed above — the resulting merger rate in 
terms of stellar mass, M, or baryonic mass, Mharvnns, can be 
enhanced near this scale ( iBundv et al.1 120091: IStewiarti l2009t 
iHopkins et al.l l2009a) and strengthen the dip feature. This 
is because the dynamical friction timescale (and hence the 
merger timescale) depen ds firstly on the mass ratio of the 
merging systems (e.g. Boylan-Kol chin et al.P 2008 and refer- 
ences therein). Low mass ratio mergers happen only after a 
long time, and one-to-one mergers happen quickly. At low 
mass, one-to-one mergers in DM correspond to minor merg- 
ers in baryons and hence the baryonic content of a halo as- 
sembles more slowly than the halo. This is due to the steep 
relation between stellar mass, Mg, and halo mass, Mh, below 
M* (see Fig. |9]l. At masses near and just above the scale 
at which Mg(Mh) flattens, more frequent minor halo merg- 
ers can host one-to-one baryonic mergers and thus the bary- 
onic assembly rates increase markedly. At even higher mass 
scales, above the flattening at M* of the Mh(Mg)-relation, the 
infalling galaxy is much more likely to become and stay a 
satellite, even if the baryonic mass ratio is close to one. At 
this point, the baryonic merger rate drops again. 

As a result, low stellar (baryonic) mass galaxies get de- 
pleted more rapidly in one-to-one mergers, and pile-up around 
the M* mass scale (M* ^ 10' 1 M©). This scenario also sug- 
gests that the two populations of galaxies might be thought 
of as central galaxies occupying the bump with an increas- 
ing fraction of satellite galaxies at lower masses that possibly 
form the second, faint, component of the mass function, lead- 
ing to a steepening of the low-mass end of the mass function. 
Such satellites in a galaxy or group halo would orbit for a long 
time before merging and hence remain visible as individual 
galaxies. 

If this enhanced assembly scenario were important, it would 
also indicate that many major mergers between star-forming 
(disk) galaxies must not lead to a final destruction of the 
disk and truncation of star formation as the bump signature 
is already apparent in the blue galaxy population at z ~ 1. 
Such mergers must produce a remnant that is still a blue star- 
forming disk galaxy. Such outcomes of mergers seem possi- 
ble if the gas fraction is large (Scrinael & Hernauist 2005; 
[Robertson elaLl 120061: iGovernato et al.1 120071: iHopkins et all 



l2009bl) . Also, we would expect the dip feature to deepen with 
time but, in contrast with SFR scenario above, remain rela- 
tively fixed with respect to mass. 

Steepening Faint-End Slope. So far we have focused on 
the dip in the stellar mass function. We now turn to the steep 
slopes at the low-mass end. To the extent that this feature is 
unrelated to the dip forming process mentioned above, we can 
consider several possibilities. To begin with, many authors 
interpret the faint-end slope as rising to the maximum set by 
the halo mass function (with a value of -2). Such behavior is 
expected if feedback processes that regulate the star formation 
efficiency become increasingly independent of Mhaio below 
some scale. At low masses, for example, supernova feedback 
may so dominate over the ability of halos to retain gas, that 
the effective gas or stellar mass fraction becomes a constant, 
independent of halo mass. 

5.3.2. Red Galaxies 

We now turn to the behavior of the red galaxy mass func- 
tion. In the redshift bins z = 0.5 and z = 0.3, where the faint 
red component of the mass function is sampled well, the num- 
ber density of all galaxies at M < 9.5 does not change much 
(Fig. [8] and Table O. Yet, the number density of faint red 
galaxies increases by 0.45 dex, a combined effect of the nor- 
malization increasing from 0f(z = 0.5) = 0.28 to 0f(z = 0.3) = 
0.49 and the characteristic mass of the low-mass component 
increasing fromM;(z = 0.5) = 9.41 to M*{z = 0.3) = 9.54. This 
evolution in faint red galaxies is nearly perfectly mirrored by 
a decrease in number of faint blue galaxies. As we have noted 
above, the faint end slope of the active and the passive popu- 
lation is consistent with being equal. 

The similarity of the steep faint-end slopes of the active 
and faint passive populations and their reciprocal change in 
number density is suggestive of the latter originating from 
the former by shutting off star formation. It is very tempt- 
ing to identify these two faint galaxy populations sharing 
the same steep faint-end slope with Sm/lm/dlrr galaxies, 
and faint spheroidal galaxies, respectively. It is established 
that passively evolving dw arf galaxies cluster around mas- 
sive galaxies (jZehavi et a l. 2005; Haines et al. 2006, 200"^ 
Carlberg et al. 2009), and tidal interactions or ram pressure 
stripping may well lead to quenching and some subsequent 
phase mixi ng turning Sm/I r n galaxies into faint spheroidal 
galax ies (iMayer e t al7'2001': 'Grebel et aLl l2003t iHaines et all 
120071) . Also, iMcCracken et al. (2008D finds fliat in the 
CFHTLS in the redshift bin 0.2 < z < 0.6 the clustering ampli- 
tude for faint red galaxies is actually higher than that of bright 
red galaxies, which is to be expected in our interpretation. 

6. SUMMARY 

Following on previous studies of the stellar mass function 
in the COSMOS field (Ilbert et al. 2009b; Pozzetti et al. 200^ 
iBolzonella et al.i 2009). we present a new analysis of this data 
set that provides mass functions to fainter limits than has been 
previously probed at z < 1. The resulting increase in dy- 
namic range allows us to characterize and study features in 
the shape of the stellar mass function that deviate from a sin- 
gle Schechter function. We have tested whether these features 
could be introduced by a variety of systematic effects includ- 
ing both catastrophic photometric redshift errors as well as 
increasing photo-z uncertainty at low masses, differences in 
the way stellar population models account for TP-AGB stars, 
and the ability of stellar mass codes to convert from luminos- 
ity to mass. We conclude that our results are robust to these 
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effects, although the data are still limited by cosmic variance. 
Our key results follow: 

• Neither the total nor the red (passive) or blue (star- 
forming) galaxy stellar mass functions can be well fit 
with a single Schechter function once the the mass 
completeness limit of the sample probes below roughly 
3 X 1O'M0. We model this more complicated behavior 
using a double Schechter function with 6 free parame- 
ters, and present the fitting results for 4 redshift bins to 
z= 1. 

• The bimodal nature of the mass function is not solely a 
result of the blue/red dichotomy. Indeed, the blue mass 
function is already bimodal at z 1 . This suggests a 
new 2-component model for galaxy formation that pre- 
dates the appearance of the red sequence. 

• We propose two interpretations for this bimodal distri- 
bution, focusing on the "dip" in the blue mass func- 
tion at ^10"^Mq. If the gas fraction increases at lower 
stellar masses, galaxies with Mbaiyon ~ 10'°Mq would 
shift to lower stellar masses, creating the observed dip. 
This would indicate a change in SFR efficiency, per- 
haps arising from the influence of supernovae feedback 
which likely sets in below scales of ^10'"Mq. In this 
picture, the baryonic mass function should not show a 
dip. Using published (cold) gas fractions as a function 
of stellar mass, we show that cold gas alone is not suf- 
ficient to eliminate the dip in the baryonic mass func- 
tion, but that the addition of the hard to detect warm gas 
could potentially flatten out the baryonic mass function 
considerably. 

• Alternatively, the dip could be created by an enhance- 
ment of the galaxy assembly rate at ^^IO'^Mq, a phe- 
nomenon that naturally arises if the baryon fraction 
peaks at Mhaio IO'^Mq. In this scenario we would 
identify the galaxies occupying the bump around 
with central galaxies and the increasing fraction of 
satellite galaxies at lower mass with the second, fainter, 
component of the mass function and in particular the 
steep faint-end slope. 



• The low-mass end of the blue and total mass functions 
exhibit a steeper slope than has been detected in previ- 
ous work. This can be interpreted as a steepening slope, 
one that may increasingly approach the halo mass func- 
tion value of -2. 



• While the dip feature is apparent in the total mass func- 
tion at all redshifts, it appears to shift from the blue to 
red population, likely as a result of transforming high- 
mass blue galaxies into red ones. 

• At the same time, we detect a drastic upturn in the 
number of low mass red galaxies. Their increase with 
time seems to reflect a decrease in the number blue sys- 
tems and so we tentatively associate them with satel- 
lite (dwarf spheroidal) galaxies that have undergone 
quenching due to environmental processes. 

While the broad dynamic range of the COSMOS data set 
allows us to begin to characterize some of the more subtle 
features of the stellar mass function, the single COSMOS field 
is still limited by cosmic variance. Future work over several 
fields will verify the trends here and shed light on possible 
interpretations by constraining how these features evolve with 
redshift. 
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